Flat-Band Ferromagnetism as a Pauli-Correlated Percolation Problem 
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We investigate the location and nature of the para-ferro transition of interacting electrons in 
dispersionless bands using the example of the Hubbard model on the Tasaki lattice. This case 
can be analyzed as a geometric site-percolation problem where different configurations appear with 
nontrivial weights. We provide a complete exact solution for the ID case and develop a numerical 
algorithm for the 2D case. In two dimensions the paramagnetic phase persists beyond the uncorre- 
lated percolation point, and the grand-canonical transition is via a first-order jump to an unsaturated 
ferromagnetic phase. 
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Introduction. — The interplay of the Coulomb interac- 
tion with the Pauli principle was already recognized by 
Heisenberg [l[ to give rise to a ferromagnetic exchange 
interaction, also encoded in Hund's rule about aligned 
spins in a partially filled shell. For a many-body sys- 
tem of correlated electrons with a flat band, when the 
interaction energy completely dominates over the kinetic 
energy, the ferromagnetic instability is one of the few 
problems for which exact results are available, albeit for 
a restricted range of fillings [l-Q • 

Flat band systems are receiving a great deal of at- 
tention right now, in particular with the view of real- 
izing new many-body phases there (see (8l-[l2| and ref- 
erences therein); in this context, the possibility of ferro- 
magnetism as a many-body instability is also being con- 
sidered It is therefore timely to provide a detailed 
study of the phase diagram and the critical properties 
of this form of magnetism: we analyze a flat-band ferro- 
magnet with an on-site Hubbard interaction of strength 
[/ > 0. For /7 = 0, any state involving electrons occupy- 
ing the flat band only is trivially a ground state. 

Crucially, this degeneracy is only partially lifted when 
a repulsive /7 > is switched on. First, since the flat 
band permits well-localized real-space electronic wave 
functions, at low density electrons can be placed on the 
lattice so that they do not overlap. Second, even if they 
do overlap, they can still avoid paying an energy penalty 
U\ the basic reason is that the Pauli principle, by de- 
manding an antisymmetric pair wave function, makes the 
overlap between two electrons on the same site vanish 
provided they are in a symmetric spin state. This is the 
origin of flat-band ferromagnetism. 

As the density of electrons increases, ferromagnetic 
clusters of increasing size appear. The degeneracy, m + 1, 
of a ferromagnetic cluster containing m electrons, gives 
differing weights to different clustering of electrons. The 



ferromagnetic transition corresponds to the emergence of 
a cluster containing a nonzero fraction of the electrons. 

An early remark by Mielke [2] likened this problem 
to one of percolation. Mielke and Tasaki [sl, Q noted 
that, for a class of flat-band ferromagnets on particular 
decorated lattices, the percolation problem in question 
is not a standard one [3, [HI but rather one including 
nontrivial weights. 

Here, we develop this analogy in detail. First of all, we 
point out that the interaction between the clusters, on 
account of its "statistical origin" in the Pauli principle, 
is unusual in that it is range-free and purely geometric — 
two particles interact only if they form part of the same 
cluster. The interaction is genuinely many-body in that 
it cannot be decomposed into a sum of pairwise terms. 
It is effectively repulsive and only depends on the size 
of the cluster, irrespective of its shape. Despite its long 
range, the statistical interaction does saturate. 

This motivates the study of the resulting unusual per- 
colation problem, which we call Pauli-correlated perco- 
lation (PCP). We find that it has a number of interest- 
ing features in its own right. It provides an instance 
of a problem in the quantum physics of strongly cor- 
related electrons which can be "reduced" to a highly 
nontrivial problem in classical statistical mechanics, on 
which an entirely different set of tools can be brought 
to bear. We first demonstrate some special features of 
this problem by providing a complete exact solution of 
the one-dimensional (ID) version of this model, which 
corresponds to a sawtooth lattice potentially realized 
in strongly correlated sawtoothlike compounds such as 
CeRh3B2 [El. Unlike standard percolation, this exhibits 
a tendency to break up large clusters as well as a de- 
velopment of spatial (anti)correlations. Its percolation 
transition at full filling is continuous. 

Next, we carry out an analysis of the phase diagram 
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FIG. 1. (Color) Left: Two-dimensional Tasaki lattice. A trapping cell contains five sites (dashed red lines). The green circles 
and lines show the ID variant of the lattice (sawtooth chain). Right: Snapshots of configurations for standard and Pauli- 
correlated percolation for small deviations from critical concentration. Panels a) and b) show snapshots (lattice extension 
£ = 200) of configurations for standard percolation for concentrations pi — 0.574 and p2 = 0.6 {pc = 0.592746...), while 
panels c), d), e), and f) show snapshots for Pauli- correlated percolation for p3 — 0.62 (paramagnetic), p4 = 0.65, ps = 0.7 
(phase-separated), and pq — 0.78 (ferromagnetic). Pink color denotes the largest cluster. 



for the two-dimensional (2D) Tasaki lattice, a decorated 
square lattice (see left panel of Fig. 1). Using a numer- 
ical algorithm custom-tailored to the problem at hand 
by extending the Hoshen-Kopelman and Newman- Ziff al- 



gorithms ll, |T3 for standard percolation, we establish 
that the ferromagnetic transition does indeed take place 
at a filling comfortably in excess of the corresponding 
well-known percolation transition on the square lattice 
at Pc = 0.592746 . . . [H, In the grand-canonical en- 
semble, this transition is of first order; in the canonical 
ensemble we find concomitant phase-separated states (see 
Fig. [U for some examples). 

Pauli- correlated percolation and flat-band 
ferromagnetism. — As a representative system with 
a dispersionless (flat) band, let us consider the Tasaki 
model [1, [1] , although our approach in principle can be 
adapted to other fiat-band lattices. The enumeration of 
all ground states of the repulsive Hubbard model on the 
Tasaki lattice maps to a percolation problem where each 
occupied site on a hypercubic lattice corresponds to an 
electron localized in a trapping cell (whose wave function 
only overlaps with that of electrons in adjacent cells.) 
(The details of this mapping, which are unimportant 
for the following, are relegated to section I Sl[) . All 
ground states can be labeled by the possible geometric 
configurations of n electrons distributed over A/" traps, 
labeled by and a nontrivial weight of each state Q 
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W{q) = l[e^\^H\a\^l) 



(1) 



which arises because of the spin degeneracy of the ferro- 
magnetic cluster of size \Ci\ in configuration q {Mq de- 
notes the number of clusters in the system). Here 
is a fugacity which can be used to tune the number of 
electrons in a grand-canonical ensemble. 

The expectation value of an operator A is given by the 
usual expression 



{A) = 



For the grand-canonical ensemble, the sum over q runs 
over all configurations of n = 0, . . . , A/* electrons while for 
the canonical ensemble, it is restricted to configurations 
with a given number of electrons n. 

From the point of view of magnetism, a particularly 
important observable is the square of the total spin S'^ 
which can be written for a particular geometric configu- 
ration q in two equivalent ways 
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(2) 



(3) 

In the first form, the contribution from each cluster is 
manifest while the second form relates it to ng(/), the 
normalized number of clusters of size /, i.e., a quantity 
which plays a central role in percolation theory [3, Ell . 

Quantum- statistical interaction. — Important differ- 
ences arise between our Pauli-correlated percolation and 
the standard one [with trivial weight factor W{q) = 1]. 
The weight factor can be cast as a pseudo-Boltzmann 
weight of statistical origin, W{q) = ex.p[\iiW{q)]. The 
resulting effective entropic interaction, \nW{q)^ has the 
following properties. First, it is repulsive — a group of 
m electrons has maximal weight 2^ if they form isolated 
one-electron "clusters," and minimal weight m + 1 if they 
form a single cluster. These extreme cases show that 
ln(m + 1) < lnW{q) < mln2 saturates, i.e., is never 
superextensive unlike other long-range interactions. Be- 
fitting its quantum statistical origin, the interaction is 
range- free — the shape of the cluster is unimportant, only 
its number of electrons matters. Note also that the in- 
teraction is a genuinely many-body one: due to the form 
\nW{q) = ln(m + 1) it cannot be written as a sum of 
two-particle terms. 

Taking all of this together demonstrates that this in- 
teraction gives rise to an entirely novel "Pauli-correlated" 
percolation problem, of interest in its relevance to flat- 
band ferromagnetism and as a physically motivated ex- 
ample of a nonstandard percolation problem with an un- 
usual weight. 

Exact solution in one dimension. — We first provide 
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FIG. 2. (Color online) Left: n{l) (top) and g{l) (bottom) at 
p = 0.99 for PCP (solid line) and standard percolation (dotted 
line) in ID. Right: Deviation of (*S'^)/<S'max from saturation 
for large p for PCP and standard percolation in 2D. 



a complete solution of the ID Tasaki model (sawtooth 
chain) [^-Q- ^ solution of the problem can be obtained 
with the help of a transfer matrix [7] despite the long- 
range nature of the statistical interaction (for technical 
details see l Sip . For a given electron density p = n/JV we 
find 



n{l) = 



4(1 -p)3 



{1^1) a' 



a = 



p 



(4) 



This cluster-size distribution (Fig. O left panel) has a 
maximum at /"^ > 1 for p > 0.8 moving along l'^ ^ —{1 -\- 
1/ In a) for p ^ 1. This is unlike the standard percolation 
result (1 — pY p^ [14], which drops monotonically with / 
and thus has a maximum at /"^ = 1. 

The macroscopic magnetic moment vanishes for p < 
1, with a continuous onset at percolation, pf = 1, as 
{pf -p)-^: 



The connected pair correlation function 



(5) 



9{\i-j\) = {niUj) 
= -(1- 
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{ni){nj) 



(6) 



yields a correlation length = —1/(2 In a) that diverges 
as {pf — p)~^ when p ^ Pf = 1. By contrast, for stan- 
dard percolation there are no nontrivial pair correlations: 
9{\i-j\) =p{l-p) Sij. 

The negative sign in Eq. (|6|) shows that the interaction 
is repulsive — electron positions anticorrelate. The long- 
range and many-body nature of the interactions leads to 
a nontrivial cluster size distribution favoring an approx- 
imately uniform spacing of vacant cells. 

The phase diagram in 2D. — The 2D case is not 
amenable to exact solution. Here we examine the 2D 
PCP numerically. Due to the nontrivial weights ([T]), sim- 
ple random sampling used for the conventional percola- 
tion is insufficient. 

SIS 



Going beyond standard numerical schemes 
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FIG. 3. (Color online) Density of electrons p vs. chemical 
potential /i, controlling the filling of the flat band in 2D. Inset: 
Histograms of density for a "finite-size" critical value /x = /ic- 



square lattices with periodic boundary conditions as fol- 
lows. In the grand-canonical ensemble, we simply choose 
a site and if it is empty (occupied), propose to insert (re- 
move) an electron. In the canonical ensemble we generate 
a new configuration q2 from the given one qi by random 
permutation of two sites in order to ensure a fixed number 
of electrons. The new configuration is accepted with the 
Metropolis probability mm[l^W{q2)/W{qi)]. In addi- 
tion, we have employed exchange Monte Carlo steps [l8| 
for the grand-canonical simulations. Clusters are identi- 
fied in two different ways: 1) using a modified Newman- 
Ziff algorithm [l^ which locally updates cluster label- 
ing for fixed number of occupied sites and 2) using the 
Hoshen-Kopelman algorithm [l6[ which makes a global 
update. 

Our central results are the following. 
The percolation transition is of first order as already 
suggested by visual inspection of individual configura- 
tions (Fig. [T]). Grand-canonical simulations exhibit a 
jump at a chemical potential fic between densities p- 
and p+, fixed by equal-sized peaks in the histograms as 
shown in Fig. [3l We estimate the jump to occur be- 
tween densities p- around 0.63(1) and p+ ~ 0.75(2). In 
between, in our canonical simulations for finite systems, 
ferromagnetism appears to set on smoothly. Figure [H 
shows {S'^)/S^g^^ [where S^^^^ = |(| + 1)] for systems 
up to 270 X 270 sites. Additionally, the cluster-size distri- 
bution n{l) indicates the emergence of a large component 
without passing through a scale-free critical distribution. 

Extent of nonpercolating (paramagnetic) phase. — The 
critical density for PCP exceeds that of the standard case 
HI, [l3), see Fig. m where the macro- 
0.62 is seen to scale to zero with 



have implemented efficient importance sampling on CxjC 



{pc = 0.592746 . 
scopic moment at p 
system size. This reflects the breakup of the large clus- 
ters due to the repulsive effective interactions. 
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FIG. 4. (Color online) Left: Finite-size scaling of the square 
of the magnetic moment (*S'^)/*S'max for the 2D Tasaki lattice 
(up to £ = 270). Right: Normalized number of clusters n{l) 
for C = 200. 



Percolating phase with unsaturated ferromagnetism. — 
For higher densities, there appears a regime of unsatu- 
rated feiiomdignetism^ where {S'^) / S^^^^ < 1, illustrated 
by p = 0.78 in Figs. [1] and HJ The existence of this regime 
is transparent from the percolation viewpoint: In stan- 
dard percolation the largest cluster excludes a nonzero 
density, N'p{l — pY for p 1, c being the coordina- 
tion number of the lattice. For PCP, with its repulsive 
interactions, this is amplified (see right panel in Fig. [2]). 
However the power law is identical to that of standard 
percolation, showing that the repulsive interactions do 
not immediately lead to a breakup of the largest clus- 
ter, presumably on account of the high entropic cost of 
arranging voids into continuous lines separating two clus- 
ters. 

We note some features of the canonical ensemble aris- 
ing in the phase-coexistence regime. The high-density 
phase appears to form first as a compact nonpercolating 
object with a macroscopic magnetic moment [the config- 
uration at j9 = 0.65 shown in Fig. [IJd) is in this region]. 
For higher densities, including the case p = 0.7, the fer- 
romagnetic phase spans across the system [compare Fig. 
[D^e)]. The details of the phase-separated regime there- 
fore contain much which is different from standard per- 
colation including the bootstrap and correlated variants 
2C 



, in particular with regard to properties which are 
of interest to fiat-band ferromagnetism. These topics are 
the subject of ongoing studies |2l|. 

Conclusions and perspectives. — We have considered 
Pauli-correlated percolation, an unusual percolation 
problem arising in a strongly correlated fiat-band sys- 
tem, where the weights of the geometrical configurations 
take nontrivial values due to the spin degeneracy and the 
Pauli principle. 

The Pauli-correlated percolation problem can be ex- 



amined exactly in ID and simulated efiiciently in 2D. 
We found that the effectively repulsive interaction leads 
to a breaking up of the clusters, and thus to a first-order 
grand-canonical transition in 2D, at a density which is 
higher than that of standard site percolation. For the 
underlying 2D Tasaki-Hubbard model our results imply 
ground-state ferromagnetism in a range of electron fill- 
ings from 0.21(1) to 1/3. 

Besides the ID realization mentioned above [5] and the 
hope of discovering corresponding quasi-2D materials, in 
[12| the possible realization of fiat-band ferromagnetism 
in organic polymers was discussed. On the other hand 
the 2D version of the Tasaki lattice is so simple that it 
seems to be a reasonable candidate for realization as a 
system of cold atoms in optical lattices lo|, 22|. 
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SUPPLEMENTARY INFORMATION 

SI. Pauli-correlated percolation and flat-band 
ferromagnetism 

These can be most transparently studied by consider- 
ing the Hubbard model, which describes electrons hop- 
ping on a lattice which interact with each other via an 
on-site repulsion. The Hubbard Hamiltonian reads 

in standard notations. 

Let us consider the Tasaki lattice [sl, Q (Fig. [ST|) as a 
representative for a 2D system with a dispersionless (flat) 
band. Tasaki's lattice decoration 0, 0| can be performed 
in arbitrary dimension (and also for other lattices). This 
allows a direct comparison of ID and 2D. For the lowest- 
energy one-electron band to be completely dispersionless, 
the two relevant hopping integrals obey the relation t' = 
\fct > 0, where c is the coordination number of the un- 
derlying lattice. The one-electron states in the flat band 
can be taken as localized on trapping cells. In 2D each 
trapping cell consists of one site of the underlying square 
lattice and four neighboring decorating sites (Fig. [ST|) . 
A localized eigenstate of energy Si = —ct = — 4t of an 
electron with spin a is given by ll^crl^), where r = x^y 
runs over the sites of the underlying square lattice, 11^^ = 
1 + c"*" 1 -^0^,1 -24^^, and 

x—^,y,a x-\-^,y,a x,y—^,a x,y-\-^,a -^ji/jcr' 

|0) denotes the vacuum state. 

Localized many-particle ground states for n > 1 elec- 
trons and U > require (i) that each cell be empty or 
singly occupied and (ii) that electrons in adjacent cells 
be in a symmetric collective spin state. As a result, the 
spin degeneracy for a contiguous cluster of m electrons 




FIG. SI. (Color online) Two-dimensional Tasaki lattice [hop- 
ping integrals t (thick lines) and (thin lines)]. A trapping 
cell contains five sites (dashed red lines). The green circles 
and lines show the ID variant of the lattice (sawtooth chain). 



is reduced from 2^ to that of a spin S = m/2 multiplet, 
m + 1. 

This is all that's needed for the mapping to PCP: all 
ground states can be labelled by the possible geomet- 
ric configurations of n electrons distributed over M cells 
(they are labelled by g = 1, . . . , ('^)), and the non-trivial 
weight of each state is given by Eq. ([1]), where Mq de- 
notes the number of separated clusters, and \Ci\ denotes 
the number of electrons in cluster i. 

We have verified by exact diagonalization that this 
set of localized many-body states spans the ground-state 
space for Tasaki lattices up to = SA/" = 48 sites. Nu- 
merical results for degeneracies are given in Table [Sl] and 
agree with the corresponding number of geometric config- 
urations. Note that even for the small lattices considered 
in this context, not all sectors up to n = A/" electrons are 
numerically accessible. In addition, the Hubbard model 
allows for higher electron fillings which are beyond the 
the geometric picture of the present paper. 



S2. Pauli-correlated percolation in one dimension 

One-dimensional PCP can be analyzed using a 
transfer-matrix method. The crucial step is to choose 
suitable representative configurations for the quantum 
states with different values of , for example by putting 
all the down spins right from the up spins in a cluster . 

We set z = exp(/i) and introduce the grand partition 
function S(z,A/') of a percolating system as the sum of 
probabilities of all possible random realizations. S(z, A/") 
will be written in terms of the transfer matrix T as 
S(z,A/') = TrT-^. The transfer matrix for the Pauli- 
correlated percolation in one dimension then reads 



T(t,0) T(t,t) T(t,;) 



\ 




1 


l\ 






z 


z 


/ 










.(8) 



The matrix elements T{ni^ni-^i) correspond to the pair 
of neighboring sites i and i -\- 1 and acquire the value 1 
{z) if the site i is empty (occupied). 
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TABLE SI. Ground- state degeneracies gjv{ri) of the 2D Tasaki model, as obtained by exact diagonalization for U = oo and n 
electrons on N = 3J\f sites. 
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10 


16 


n 


1 2 3 4 5 


12 3 4 


12 3 4 


En/t 


-4 -8 -12 -16 -20 


-4 -8 -12 -16 


-4 -8 -12 -16 


9M{n) 


16 96 256 372 336 


20 160 640 1380 


32 448 3584 18008 



In order to determine the results in terms of p rather 
than we first calculate the average occupation number 
of the site (n^) which should be equal to p. Thus we have 



{Ui) 



TrT-^N 
TrT-^ 



■P. N: 




(9) 



In what follows we consider the thermodynamic limit 
M ^ OO. Using the left and right eigenvectors of T, one 
obtains p{z) = 1 + 4 ^ — ^/l -\- / [1 + 4 2;] , and inversion 
of this function leads to z{p) = p{2 — p)/[4{l — p)^]. 

We turn to the calculation of the average number of 
clusters of size / (normalized by the lattice size A/") n{l) 
[1^ . To fix the cluster of length / we start with an empty 
site, then we have a string (cluster) of / occupied sites, 
and the last site of this string is followed by an empty 
one. To calculate n{l) we have to replace the product of 
a sequence of / + 1 T- matrices by the product SC^~^F, 
where 





(0 


1 


1^ 




(0 





0) 









o\ 


s = 
















z 


z 




z 















0) 













u 





0/ 



That yields 



n(l) 



TrT- 



sc 



TrT-^ 

After straightforward calculations we arrive at 



(10) 



(11) 



p 



(12) 



(2-p)2 ^" ' ' ' ' 2-p 
see Eq. (Hf. In these calculations we have used the rela- 



tion 




(13) 



Next, we use the transfer-matrix approach to calcu- 
late the (pair) site-occupation correlation function g{l) = 
{niUi^i) - {ni){ni^i) = {uirii^i) - p^ . Using the matrix 
N defined in Eq. ^ and calculating TrT-^-^NT^N we 
get 

g{l) = -(l-pfa^\'\. (14) 
Finally, we calculate the pair connectivity r(n,n + /) 
(the probability that two sites n and n + / are both oc- 
cupied and belong to the same cluster). For this purpose 
we have to consider the quantity 



r(n,n + /) 



Tr(T-^-^NCO 



p- 



TrT-^ 



2-p 



(15) 



For completeness we mention that the transfer-matrix 
approach can also be applied to the standard percolation 
in one dimension. In this case we have 



S = 




C = 




(16) 



and formulas ([9j), ([TT]), and ([l5]) yield z = p/{l 
n(l) = (1 -pfp\ g(l) =p{l -p)Si^o, and r(n,n + /) = 
p - p\ respectively. These results can of course also be 
obtained by more elementary means and are well-known 
in standard percolation theory Q . 

In Fig. [2] we illustrate n{l) and g{l) for both types of 
percolation at p = 0.99. 



